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ABSTRACT 

Aims. We have analysed optical spectra of BL Lacertae, the prototype of its blazar subclass, to verify the broad Ha emission line 
detected more than a decade ago and its possible flux variation. We used the spectroscopic information to investigate the question of 
the BL Lacertae parent population. 

Methods. Low- and high-resolution optical spectra of BL Lacertae were acquired with the DOLORES spectrograph at the 3.58 meter 
Telescopio Nazionale Galileo (TNG) during four nights in 2007-2008, when the source was in a relatively faint state. In three cases 
we were able to fit the complex Ha spectral range with multiple line components and to measure both the broad Ha and several 
narrow emission line fluxes. 

Results. A critical comparison with previous results suggests that the broad Ha flux has increased by about 50% in ten years. This 
might be due to an addition of gas in the broad line region (BLR), or to a strengthening of the disc luminosity, but such flux changes 
are not unusual in Broad Lined active nuclei. We estimated the BL Lacertae black hole mass by means of its relation with the bulge 
luminosity, finding 4-6 x 10* Mq. The virial mass estimated from the spectroscopic data gives instead a value 20-30 times lower 
An analysis of the disc and BLR properties in different AGNs suggests that this discrepancy is due to an underluminosity of the BL 
Lacertae BLR. Finally, we addressed the problem of the BL Lacertae parent population, comparing its isotropic quantities with those 
of other AGN classes. From the point of view of the narrow emission line spectrum, the source is located close to low-excitation radio 
galaxies. When one also considers its diffuse radio power, an association with FR I radio galaxies is severely questioned due to the 
lower radio luminosity (at a given line luminosity) of BL Lacertae. The narrow line and radio luminosities of BL Lacertae instead 
match those of a sample of miniature radio galaxies, which however do not show a BLR. Yet, if existing, "misaligned BL Lacertae" 
objects should have entered that sample. We also rule out the possibility that they were excluded because of a QSO optical appearance. 
Conclusions. The observational constraints suggest that BL Lacertae is caught in a short term transient stage, which does not leave a 
detectable evolutionary "trace" in the AGN population. We present a scenario that can account for the observed properties. 
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1. Introduction 

According to the commonly accepted scenario, the central en- 
gine of active galactic nuclei (AGNs) is a supermassive black 
hole (SMBH) fed by infall of matter from an accretion disc. 
Inner fast-moving clouds produce broad emission lines, which 
may be obscured by absorbing material, while outer clouds 
are responsible for n arrow emission lines. A bout one fifth of 
AGNs is radio-loud dKeUermann et alj|1994l) . showing plasma 
jets sometimes extending on Mpc scales. Among them, BL Lac 
objects and flat spectrum radio quasars (FSRQs) form the blazar 
class, characterized by variable emission from the radio to the y- 
ray band, with flux variations on time scales from hours to years, 
high radio and optical polarization, core-dominated radio mor- 
phology, flat radio spectra, and apparent superluminal motion 
of radio jet components. Their properties are explained in terms 
of plasma relativistic motion in a jet pointing at a small angle 
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with the Une of sight, with consequ ent beaming of the observed 
radiation dBlandford & Reeslll978h . Hence, the continuum radi- 
ation of blazars is dominated by the relativistically beamed non- 
thermal radiation from the jet. 

In FSRQs, thermal emission from the disc may be observable 
in the optical-ultraviolet band when the source is not in a flaring 
state. Disc signatures were d etecte d e.g. in 3C 273 (j Smith et al] 
1993'; von Monti gnv et all '1997; Grandi & Palumbo ,2004^ 
Turleretal. 2006), 3C 2 79 (Pianetal. 1999), PKS 1510- 
089 (Kataokaetal]_p^ [commando etal. 2OO 90, 3C 34 5 
(Bregman et aL 198 0). and 3C 454.3 (Raiteri et al..2007bll2008l) . 
Moreover, strong broad and narrow emission lines are usually 
present in their spectra. 

As for BL Lac objects, according to the original defini- 
tion, they may show at most weak emission lines, with equiv- 
alent widths not exceeding 5 A in the rest frame (iStickel et^D 
1991). This seems to be due not so much to low line fluxes , 
but rather to a high continuum flux jScarpa & Falomolll997h . 
Indeed, strong emission lines, in particular broad ones, have oc- 
casionally been detected in the spectra of BL Lac objects in faint 
states. These include BL Lacertae, the prototype of the blazar 
subclass named after it (IVermeulen et al.l I1995L ICorbett et al.l 
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1996, 



1987 



2000h. and the distant source AO 0235+ 164 (ICohen et alJ 



Nilsson et al ]|l996HRaiteri et al.ll2007ah 



The unified scheme for radio-loud AGNs predicts that BL 
Lac objects and FSRQs are the beamed counterparts of Fanaroff- 
Riley type 1 (PR 1) and Fanaroff-Riley type 11 (PR 11) radio galax- 
ies, respectively, even if these corresp ondences were qu estioned 
by various observing evidences (e.g. iTadhunte^ l2008l) . In par- 
ticular, many BL Lac objects show hig h, PR 11-hke extended 
radio powers and morphologies (see e.g. lLandt & Bignailll2008l 
and references therein). The distinction between blazars and ra- 
dio galaxies, as well as that between PSRQs and BL Lac ob- 
jects, has been discussed by several authors, and different cri- 
teria h ave been proposed , involving the value of the Ca H&K 
break (iMarcha etalJll996h or the strength of the oxygen-narrow 
emission-lines dLandt et alj|2004l) . 

A powerful way to classify AGNs is their position in 
diagnostic diagr a ms comparing selec ted emission line ratios 
(lHeckmanl[l980t iBaldwinet aLfll98lh . In particular, ratios of 
lines close in wavelength, Uke [0111]/HyS, [Nll]/16583/Ha, 
[S11]AA6716,6731/Ha', and [Ol]/Ua are expected to be the 
most reliable ones ( Veilleux & Osterbrock 1987). The applica^ 
tion of diagnostic diagrams to radio-loud galaxies bv iLaing et al.l 
11994) confirmed former suggestions that PR 11 sources can 
be divided between high-excitation g alaxies (HEG) and low- 
excitation galaxies (LEG). In particular.'Buttigl ione et al.l (1201 Ot) 
verified this dichotomy when analysing the radio sources be- 
longing to the well known 3CR catalogue. They found promi- 
nent broad lines in a sub-sample of HEG, but not in LEG. 
Moreover, they saw that HEG are associated with very power- 
ful PR 11 only, while LEG are spread on a wide range of radio 
powers, and can be of both PR 11 and PR 1 type. Actually, the 
situation is even more complex, as the existence of miniature 
radio galaxies, characterized by extremely low radio power, rel- 
atively luminous narrow em ission lines, and no BLR, demon- 
strates (Baldi & Capetti ''2Q09). 

An analysis of the spectroscopic properties of blazars to un- 
derstand their relationship with the radio galaxies (and other 
AGN classes) is not an easy task, because the dramatic vari- 
ability of the non-thermal continuum flux strongly affects the 
appearance of lines, especially in BL Lac objects. However, this 
analysis can help clarify whether blazars differ from radio galax- 
ies only for their orientation with respect to the line of sight, or 
if they are intrinsically different sources. 

We present spectroscopic observations of BL Lacertae car- 
ried out in 2007-2008 with the 3.58 m Telescopio Nazionale 
Galileo (TNG) on the Canary Islands, to address the problem 
of its parent population. In the same period BL Lacertae was 
the target of a multiwavelength campaign by the Whole Earth 
Blazar Telescop43 (WEBT), also involving three pointings by 
the XMM-Newton sate llite. The results of the WEBT campaign 
have been presented bv lRaiteri et al.l (I2OO90 . The source was ob- 
served in a relatively faint state at all wavelengths, and a UV 
excess was clearly visible in the source spectral energy distri- 
bution (SED), which was interpreted as the signature of thermal 
radiation from the accretion disc. 



2. Spectroscopic observations and data reduction 

All spectra were taken with the Telescopio Nazionale Galileo 
(TNG), a 3.58 m optical/infrared telescope located on the Roque 
de los Muchachos in La Palma Canary Island (Spain). The ob- 
servations were made with the DOLORES (Device Optimized 
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Fig. 1. The behaviour of the Ho- broad emission line flux as a 
function of time (top panels) compared to the continuum flux 
evolution (bottom panels, /?-band magnitudes) in 1995-1997 
(left) a nd 2007-2008 (right). Blue filled squares display the re- 
sults of ICorbett et al.l (HOOO); blue empty squares represent data 
with revised flux calibration; blue diamonds show the results ob- 
tained in this paper by analysing the TNG data. Grey dots rep- 
resent data fro m the WEBT collaboration dVillata et al.l l2004t 
iRaiteri et al.l20 09). The solid horizontal lines in the upper panels 
indicate the average Ha flux in the corresponding period, while 
dotted lines are drawn through one standard deviation from the 
mean value. 



for the LOw RESolution) spectrograph installed at the Nasmyth 
B focus of the telescope. The detector was a 2048 x 2048 pixels 
E2V 4240 back-illuminated CCD, with a pixel size of 13.5;um, 
and a scale of 0.252" pixel ', which implies a field of view of 
8.6' X 8.6'. 

The spectroscopic observations were performed in service 
mode on 2007 July 10-11, September 7, and October 10, 
and on 2008 January 8, during a multiwavelength campaign 
by the WEBT; in particular, the first and last observations 
were contemporane ous to XMM-Newton pointings at the source 
(IRaiteri et al.ll2009t) . Pigure [T] (bottom right) shows the /?-band 
light curve obtained by the 30 optical telescopes of the WEBT: 
the source showed a noticeable variability, with short-term 
(intra-day or inter-day) flickering superposed on a longer-term 
trend. Blue diamonds display the results of aperture photom- 
etry performed on TNG images acquired with the Cousins' R 
filter just before/after the spectra, with the same prescriptions 
used to derive the WEBT magnitudes, i.e. an aperture radius of 
8", and background derived fro m a surrounding an nulus with 
10" and 16" radii. In this way, IVillata et aP (l2002l) estimated 
that the measured flux density includes 60% of the host galaxy 
flux. Table [U reports the /?-band magnitude^ derived from the 
TNG imaging frames according to the above WEBT prescrip- 
tions. Pigure [U shows that the photometry on the TNG images 
agrees excellently with the WEBT data and helps to put the spec- 
troscopic observations into context. 

These were performed with a 2"-width long-slit, which was 
aligned along the parallactic angle to minimize light losses due 
to atmospheric dispersion. Por each observing epoch, we ob- 
tained both low-resolution spectra with the LR-B grism, cover- 
ing the wavelength range 3500-7700 A with dispersion 2.52 A 
pixels"' and resolution 1200 (for a 2" slit), and high-resolution 



' jhttp : //www ■ oato . inaf / it/blazar s/webt/ j 



^ Throughout the paper we report uncertainties at Icr confidence 
level. 
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Fig. 2. Left panels: low-resolution TNG spectra. The position of the emission lines is indicated as well as that of the telluric 
absorption lines. Blue horizontal lines in the bottom panel mark the wavelength intervals used to calculate the D„(4000) spectral 
index defining the strength of the Ca H&K break. Right panels: high-resolution TNG spectra both before (black) and after (red) 
correction for the telluric absorption lines around the Ha line. The corrected spectra are shifted in flux for sake of clarity. 



spectra with the VHR-R grism, sensitive to wavelengths from 
6100 to 7700 A with dispersion 0.70 A pixels"' and resolution 
5000. For each grism, two consecutive frames were taken and 
subsequently combined, for a total exposure of 600 s in LR-B, 
and 1400 s in VHR-R. The resulting signal-to-noise ratios of the 
co-added spectra are ~ 50-100 and ~ 100-200 for the low- and 
high-resolution spectra respectively. 

The use of both LR-B and VHR-R grisms allowed us 
to cover the spectral range where the most relevant emission 
lines of the optical spectrum are expected, in particular the 
key diagnostic Hnes H;6, [O III]^^4959,5007, [O I]^^6300,64, 
Ha, [N II]^^6548,84, and [S II]^^6716,31. The high-resolution 
spectra were needed to disentangle Ha from the [N II] doublet. 

The longslit package of the NOAO's Image Reduction 
and Analysis Facility (IRAFj3 was used to perform bias- 
subtraction, flat-field correction, wavelength calibration, back- 
ground subtraction, spectra extraction (performed over 2" in the 
spatial direction), and flux calibration. For the last task we used 
the spectrophotometric standard star BD -i-28°4211 (Oke 1990). 
We discarded the October spectra because of their low quality. 

The absolute flux calibration of the spectra was checked 
through synthetic aperture photometry on the images taken in the 
Cousins' R band just before/after the spectra. By setting an aper- 



^ IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 



ture size of 2", we derived the flux densities at A^jf - 64 10 A 
reported in Table [1] The comparison with the spectra continuum 
around 6400 A revealed that all spectra flux densities were in 
fair agreement with the fluxes inferred from photometry, with 
the only exception of the VHR-R spectra of July 2007. Indeed, 
the scaling factors to apply to the LR-B spectra were: 0.96, 0.95, 
1.03 for July 2007, September 2007, and January 2008, respec- 
tively, while for the VHR-R spectra we obtained: 1.27, 0.99, and 
1.02. For the LR-B spectra, which extend over almost all the 
wavelength range covered by the broad-band R filter, these cal- 
ibration factors were checked by convolving the BL Lacertae 
spectra with the transmission curve of the TNG 7? -band filter, 
multiplied by the CCD quantum efficiency curve. We verified 
that the uncertainty on the calibration factors is of the order of 
2%. 

By following'Scarpa et alj (l2000l) . who gave the host galaxy 
brightness R - 15.55 + 0.02 and its effective radius re - 
4.8", and adopting a De Vaucouleur brightness profile, we es- 
timated that within 2" the host galaxy contribution is ~ 0.32 x 
10"'^ erg cm"^ s"' A"', i.e. ~ 10% of the observed flux in 
July, and ~ 8% and ~ 7% in September and January, respec- 
tively. This result is confirmed by the analysis of the Ca H&K 
break strength in the spectra. In the January spectrum we found 
that the ratio between the fluxes in the rest-frame wavelength 
ranges 4000-4100 A and 3850-3950 A is D„(4000) ~ 1.06 
(Baloghetal. 199 9i). Assumin g D„(4000) ~ 2 for the host 
galaxy (.Kauffmann et al.l l2003h . as typical of giant early-type 
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Fig. 3. Top panels: enlargement of the high-resolution, telluric- 
corrected spectra of Fig. |2] showing the results of the spectral 
fitting: red continuous lines indicate the fit to the broad and nar- 
row components of the Ha line; blue dotted-dashed lines show 
the fit to the [N 11] doublet; the green line displays the overall 
fit, which includes a tentative fit to the [S 11] line. Bottom pan- 
els: residuals of the spectral fit. Fluxes and residuals are given in 
10"'^ erg cm"^ s"' A"' units. 



Table 1. Results of aperture photometry on the TNG imaging 
frames. 



Date 






^« 




[maj 




[lO^'^'ergcm-- s"' A"'] 


2007 Jul 10-11 


14.46 + 


0.01 


3.16 + 0.01 


2007 Sep 7 


14.21 + 


0.02 


4.12 ±0.09 


2007 Oct 10 


13.79 ± 


0.01 


6.31 ±0.05 


2008 Jan 8 


14.13 + 


0.03 


4.42 ±0.11 



" i?-band magnitudes obtained from photometry with 8" aperture 
radius 

* Flux densities at 6410 A derived from photometry with 2" aperture 
radius, used to normalise the spectra. 



galaxies, this implies a host galaxy contribution not greater than 
a few percent. Thus we did not correct the spectra for the host 
galaxy contribution, because it is negligible for our purposes. 

The VHR-R spectra were corrected for the telluric absorp- 
tion bands, a particularly important step, because the 6870 A 
oxygen B-band affects the blue wing of the Ha line. This was 
done by dividing the spectra by a template obtained from the 
spectrophotometric standard star. 

Figure |2] displays the low-resolution (left panels) as well as 
the high-resolution (right panels) spectra, the latter both before 
and after correction for atmospheric absorption around the Ha 
line. 

We measured the emissi on line inten sities by means of the 
specfit package in IRAF (lKrissl[T994l) . To reduce the num- 
ber of free parameters we fixed the relative wavelength dis- 
tance between lines and required the FWHM to be the same 
for all the narrow lines. From the low-resolution spectra we de- 
rived the [O 111]/1/14959,5007 flux, and an upper limit to the Hy6 
flux. In the high-resolution spectra, we fitted the 6700-7350 A 
spectral region with seven components: a Gaussian profile for 
each of the lines: Ha (broad and narrow), [N 11]/1/16548,84, and 
[S ll]/l/16716,31,plus a linear component for the continuum. The 
results of this spectral fitting are shown in Fig. |3] The fit to the 
[S 11] /i/l67 16,6731 doublet must be considered with caution, be- 
cause this line is strongly affected by atmospheric absorption 
and thus its measurement strongly depends on the telluric cor- 
rection (see Fig.|2]i. The FWHM of the narrow lines ranges from 
~ 500 to ~ 600 km s ', while that of the broad Ho- component 
is 4200-5000 kms Including a further component to fit the 
[O I]/l/l6300,64 line resulted in a marginal detection only. 

Table|2]reports the line fluxes after correction for the Galactic 
extinction, adopting Afl =1.42, E{B — y)=0.35, and the extinc- 
tion law of lCardelli etaP (Il989h . 

3. Discussion 

3.1. Variability of the broad Ha emission line 

The first detection of broad emission line s in the optical 
spectrum of BL Lacertae was reported by IVermeulen et alJ 
(1995). They took two spectra in May and June 1995, and 
measured Ha with FWHM of 3400-3700 kms ' and flux of 
4.4 X 10"''* erg cm"^ s"', while for H/3 they obtained FWHM = 
4400km s"' and F - 1.3 x 10"'"*ergcm"^s"'. They stressed 
that such a luminous Ha line should have been observed before, 
while it is not recognizable in earlier spe ctra taken in 1975-197 6 
by Miller & Hawlev ( 1977), in 1985 by [Lawrence et all (Il996l) . 
and in 1989 by Sticke l et all (Il993l) . They estimated that the 
broad Ha flux may have increased by at least a factor 5 since 
1989. 

Soon after, ICorbett et alJ d 19961) analysed two other opti- 
cal spectra acquired in June 1995, and confirmed the results of 
IVermeulen et alJ (Il995h . They discussed the appearance of the 
Ha line as due to an increase of either the amount of gas in 
the broad line region (BLR) or the strength of the photoionis- 
ing source. In the latter case, the accretion disc seemed the most 
likely source of photoionising radiation. The signature of the ac- 
cretion disc was later recognized bv lRaiteri et alJ (l2009h as a UV 
excess in the broad-band SEDs of BL Lacertae built with con- 
temporaneous low-energy data from the WEBT and high-energy 
data from the XMM-Newton satellite. 

The spectroscopi c monitoring of BL Lacertae was contin- 
ued by ICorbett eTalJ (|2000), with eight spectra taken between 
June 1995 and December 1997, in seven of which they were 
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Table 2. Emission lines fluxes after correction for Galactic extinction. All fluxes are 10 erg cm"^ s ' ; FWHM are km s"' . 



Obs 


H/3 


[0 1111^4959,5007 


HQ'(narrow) 


[N 111^6548,84 


[S 111^^6716,6731 


HQ'(broad) 


FWHM 


2007 Jul 10-11 


<7.6 


4.7 ± 1.2 


2.4 ± 0.5 


5.3 ±0.6 


<3.1 


33.5 ±0.9 


4600 


2007 Sep 7 


< 3.0 


3.4 + 0.5 


4.8 ± 1.5 


6.1 ± 1.0 


3.9 + 0.9 


40.0 ±2.1 


5000 


2008 Jan 8 


< 2.0 


2.5 ± 0.4 


2.5 ± 0.6 


4.7 ± 0.7 


4.4 + 0.8 


34.5 ± 0.9 


4200 



able to measure the Ha broad emission line. Their results are 
reported in Fig. [1] where the upper left panel shows their line 
flux (corrected for extinction and host galaxy contamination) as 
a function of time, while the lower left panel displays their con- 
tinuum estimates superposed to the light curve by the WEBT 
collaboration ( Villata et al. 2004). They could not distinguish 
the Ho- line in June 1997, when the source was experiencing 
a big optical outburst. When considering their estimates of the 
Ha broad emission line flux together with ours, we find an av- 
erage flux of 30.6 X 10"'^ erg cm"^ s"', and a standard deviation 
of 8.8 X 10^'^ergcm^^s"'. This would imply that if we take 
into account the errors, all Ha fluxes are within one standard 
deviation from the mean value, wit h the only exception of the 
December 1997 value. Indeed, Corb ett et al.l (2000) concluded 
that in 1995-1997 there was no "compelling evidence" of "any 
significant variation", and now we could add that after about 10 
years the Ha flux has still roughly the same intensity. 

However, the WEBT photometry allows us to improve the 

ICorb ett et al. (2000) absolute flux calibration. The spectra of 
August 24, 1995, November 14, 1997, and December 7, 1997 
are the ones for which we both have cont emporaneous WEBT 
data and found significant deviations of the lCorbett et al.l (l2000h 
measures from the WEBT ones. The blue empty squares in Fig. 
[Urepresent data rescaled; the flux rescaling factors are 0.6, 0.55, 
and 0.7 for the three spectra, respectively. Thus, if we consider 
all the Hff broad emission line fluxes in the 1995-1997 period, 
rescaled to match photometric values, we find an average value 
of (23.0 + 6.6) X lO"'*" erg cm"^ s"' . This can be compared to the 
average value of (36.0 + 3.5) x 10"'^ erg cm"^ s"' obtained by 
considering our three measurements in 2007-2008, suggesting 
that a ~ 50% increase of the broad Ha intensity may have oc- 
curred in a roughly 10 year time interval. However, this is not 
unusual and it is not necessarily due to an evolutionary trend. 
Indeed, oscillations of the Ha flux u p to ~ 77% on a tim e scale 
of a few years have been reported bv lKaspi et al.l (l2000l) for PG 
quasars. 

Due to its greater distance from the central engine, the nar- 
row line region (NLR) is expected to react to variations of the 
disc luminosity on muc h longer timescales than those charac- 
terizing the BLR. Indeed ICorbett et al.l (Il996h reported fluxes of 
the narrow Ha line component and of [N II] that are inside the 
range of our results (see T a ble ID . They also identified [O I] and 
[Fe VII] . IVermeulen et al.l (Il995h measured the narrow [O III] 
and [Nil] emission lines, with fluxes simila r to ou rs. [OIII] was 
pr eviously measured also by [Sdckel et al.' ( Il993h in 1989, and 
by 'Lawr ence et al.l (1199 6) in 1985, who estimated a flux within 
a factor of 2 with respect to our spectra. In conclusion the NLR 
luminosity did not undergo significant changes in the last twenty 
years, although it must be noted that the narrow line fluxes are 
known with larger uncertainties than the broad Ha. 

3.2. Black hole mass, BLR, and accretion disc properties 

The masses of SMBH powering AGNs can be estimated in a 
number of different ways. 



One method relies on the presence of a relation between 
SMBH mass and near-i nfrared bulge luminosity as derived by 
iMarconi & Hunt! (l2003l) . Fr om the BL L acertae host galaxy 
brightness R = 15.55 + 0.02 dScarpa et al.l ll2000). correcting for 
Galactic extinction and using a mean colour index for elliptical 
galaxies of R - K - 2.71 (Mannucci et al. 2001) we obtain an 
apparent magnitude K = 11.95. This translates into an absolute 
magnitude M/c - -25.33, a nd hence \Qg{LKl Lq - 11.44. 
According to lMarconi & Hu nt (2003), this value implies Mbh ~ 
6 X IO'^Mq. 

A revision of the SMBH mass v ersus bulge luminosi ty rela- 
tion for AGNs has been proposed bv lBentz et alJ (l2009bh . using 
a sample of 26 objects observed by the Hubble Space Telescope, 
and for which SMBH masses have been estimated through the 
reverberation-mapping technique. By adopting for BL Lacertae 
Lv = 5.55 X 10'"Lo we derive Mbh = 3.76+^;^ x lO^M©, in 
substantial agreement with the previous estimate! 

Another possible approach is to use the spectroscopic infor- 
mation, and to calculate the virial mass: Mbh = /^BLRO"ji„e/G, 
where / is a factor depending on the BLR structure, kinemat- 
ics, and orientation; 7?blr is the size of the BLR; crune is the 
line dispersion , and G is the gravitational constant. Following 
Pet erson et all (l200l . / = 5.5, and cr,i„e = FWHM/2.355. 
Lacking a measurement from reverberation mapping, the size 
of the BLR can be derived from the scaling re lationship with its 
luminosity as discussed bv'Kaspi et al.l(l2005h . According to our 
data, the luminosity of the broad Ha line is ~ 4 x 10"" erg s"'; 
adopting a flux ratio H a/H/? ~ 3 (an a ssumption supported 
by the measurements by I Vermeulen et al.i |1995.) , this leads to 
/^BLR ~ 5 It day. Taking FWHM = 4600 km s"' from Table |2] 
we obtain Mr h ~ 2 x lO' Mq. Following the prescription of 
Marconi et al. (l2008h to account for the possible role of radia- 
tion pressure on the BLR, we derive a relatively small correction 
term to the black hole mass of ~ 0.6 x 10'' M©. The virial es- 
timate of Mbh is then about 20-30 times lower than the values 
derived above. 

At this stage the reason for the discrepancy is unclear. 

Further clues on the properties of BL Lacertae come from the 
relationship of the observed BLR with the accretion disc. The 
Optical Moni tor XMM-Newton observations of BL Lacertae 
presented in Raiteri et al.l (l2009l) revealed a sharp up-turn in the 
SED ultraviolet region, the characteristic signature of the Big 
Blue Bump associated with the emi s sion o f the accretion disc. 
^From Figs 6 and 7 in iRaiteri et~aL I (l2009l) one can estimate a 
disc luminosity at 2500 A of log Lv(2500A) ~ 29.4 erg s"' Hz"'. 
This value is similar to that measured in type 1 AGN (both 
radio-loud and radio-quiet) Hke Fairall 9, PG 1229+204, 3C 
390. 3, 3C 120, Mkn 509, PG 213 0+099, PG 0844+349, and Akn 
120 dVasudevan & Fabianll2009l) . The luminosity of the broad 
HyS lines of these sources ranges from 0.2 to 0.6 x lO^^* erg s ' 

The SMBH m ass of BL Lacertae was also calculated by 
'Wo o & Urnl ( l2002h through the correlation between the SMBH mass 
and the stellar velocity dispersion cr; the estimate of a is inferred from 
the m orphological paramete rs of the host galaxy, via the fundamental 
plane jJorgensen et al J 1996h . They found Mbh = 1-70 x 10** M©. 
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(iKaspi et al.ll2005h . For BL Lacertae, considering its broad Ho- 
line luminosity and the Ha/H/3 flux ratio (see above), we obtain 
a factor of ~ 15^0 lower. 

A compariso n with broad- line radio galaxies in the 3CR 
sample iButtighone et al.ll2010h gives a similar result: the ratio 
between the broad Ha and the rest-frame UV flux ranges be- 
tween 100 and 250 A, while this ratio for BL Lacertae is ~ 0.3A. 
This difference is preserved considering an AGN of very low 
BLR luminosity, ~ 100 smaller than BL Lacertae, the LINER 
galaxy NGC 4579 jBarth et alJIl 99611200 Ih . for which this value 
is ~ 275A. Apparently, the BLR of BL Lacertae is strongly un- 
derluminous with respect to its disc emission when compared to 
other AGN. 

This suggests a possible interpretation for the different val- 
ues of SMBH mass found above. In fact, Mbh scales with the 
BLR luminosity as Mbh ^blr'^- account for an underesti- 
mate of the SMBH mass by a factor of 20-30, the BLR should 
be underluminous by a factor of 70-300, in broad agreement 
to what we derived considering the ratio between BLR and UV 
fluxes. 

3.3. The parent population 

We here explore how BL Lacertae would look like when seen 
with its jet pointing at a larger angle from our line of sight, 
i.e. what extragalactic sources might represent the misoriented 
parent population. We thus consider all quantities that are not 
affected by beaming and that might eventually depend on ori- 
entation only if there is selective obscuration, as e.g. due to a 
flattened circumnuclear dust structure, e.g. a torus. 

The host of BL Lacertae is a giant elliptical galaxy of abso- 
lute magnitude Mk = -25.33 (see Sect l3.2l i. nearly two mag- 
nitudes brig hter than the characteristic absolute magnitude M* 
in this band ( Schechte^ll976t iHuang et al.ll2003l) . The emission 
from the accretion disc, to which we associate the emission in the 
near UV band at 2500 A, amounts to 10^^ "* erg s"' Hz"' . In order 
to estimate th e contrast between the disc and the host we used the 
iRaiteri et al.l 12009) black-body fit to the accretion disc emission, 
which gives an /?-band flux of 0.95 x 10 e rgcm~^ s~' A~' Th e 
host-galaxy observed magnitude = 15.55 (IScarpa et al.ll2000h . 
after correcting for Galactic extinction, translates into a flux of 
3.07 X 10"'^ ergcm"^ s"' A"'. This implies a ratio between host 
and nucleus of ~ 3, simil ar to that measured in Seyfert 1 and 
Broad Line radio galaxies dBentz et al.ll2009al) . 

The extended radio emission is also unaffected by orienta- 
tion. B y means of VLA maps at 20 cm, Antonucci & Ulvestad 
dun" estimated an extended radio flux of 40mJy. More re- 
cent VLA observa tions at 20 cm for the MOJAVE project 
(ICooper et al.ll2007l) resulted in an extended flux of 18 mJy, es- 
sentially confirming the earlier results. These high spatial res- 
olution observations might, in principle, have missed diffuse 
low-surface brightness, steep-spectrum emission. We then con- 
sider the 74 MHz flux density measured for BL Lacertae in the 
VLA Low-frequency Sky Survey (VLSS. ICohen et al.l2007h . i.e. 
Fj4 - 1.46 Jy. Even in the assumption that the radio core does 
not contribute significantly at this low frequency, this translates 
into an upper limit of ~ ISOmJy at 1400 MHz (having adopted 
a radio spectral index of 0.7). This is an extremely conservative 
limit, c onsidering that its r adio core has a typical flux of ~ 2 Jy 
(see e.g.'Kharb et al."2010'). Thus we are confident that a value of 
40 mJy at 1400 MHz is well representative of the total extended 
emission in BL Lacertae. This corresponds to a radio luminosity 
power of log Pext = 30.57 erg s"' Hz"'. 
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Fig. 4. The diagnostic index log[0 III]/[N II] vs. 
log[OI]/[NI I] of BL La certae derived from our TNG data 
and those of ICorbett et al.l ([1996) (black dot). The red squares 
and blue circles show the distributio n of LEG and HEG, re - 
spectively, according to the results bv IButtighone et all (l2010l) . 
The crossed circles mark HEG that also show a BLR. Contours 
represent the densi ty of SDSS emission line galaxies from 
iKewlev et al.](l2006h . 



/,From the point of view of the emission lines, we derived 
a broad Ho- luminosity of ~ 4 x 10"*' ergs"' (see Sect l3.2b . 
Considering the narrow emission lines, from Table |2] we esti- 
mate a [O III] luminosity of ~ 4 x lO'^'^ergs"'. The narrow 
emission line ratios can contribute to characterize the properties 
of an AGN, but unfortunately accurate measurements from our 
data are available only for [O III] and [N Il]f[ As discussed in 
Sect. 3.1, the narrow emission lines do not seem to have var- 
ied significantly in the last 20 years, so we can complement 
ou r data with those from t he literature. However, the spectrum 
of IVermeulen et al.l (Il995h does not improve the situation, be- 
cause the HjS line is detected, but a decomposition into nar- 
row and br oad component coul d not be performed. The obser- 
vations by ICorbett et al.l (Il996h cover only the red part of the 
spectrum (A > 6000A) where they saw a rather well-defined 
[O I]/16300 line. This enables us to locate BL Lacertae in a non- 
standard diagnostic plane defined by the ratios [O I]/[N II] and 
[O III]/[N II] (Fig. |4]i. As a comparison, we show in this dia- 
gram the location of the 3CR sou rces (limited to a redshift of 
0.3) from lButtiglione et alj (l2010 l). The two main spectroscopic 
classes of Low and High Excitation Galaxies (LEG and HEG, re- 
spectively) define two separate sequences, which are also present 
in the SDSS emission line galaxies from Kewley et al. (200^. 
BL Lacertae falls in a region not well populated by 3CR sources, 
but it is closer to LEG than to HEG; furthermore it lies on the 
branch of the LINERs from the SDSS. This suggests a tentative 
identification as a Low Excitation Galaxy from the point of view 
of its narrow emission line spectrum. 

Which class of objects share the properties of BL Lacertae 
described above? Let us start considering the radio-galaxies in 
the 3CR sample, see Fig. |5] The extended radio luminosity of 
BL Lacertae is at the faint end of the FR I in the 3CR, with only 



' The narrow component of Ha is heavily blended with the brighter 
broad emission; the [S II] doublet is significantly uncertain as it falls 
deep into a telluric band; our data provide us with an upper limit on Hj8 
and only with a tentative detection of [O l]/i6300. 
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Fig. 5. [O III] luminosity [erg s"'] as a function of radio lumi- 
nosity at 1400 MHz [erg s ' Hz ']. The large black dot repre- 
sents the location of BL Lacertae. Blue circles are HEG (crossed 
circles are broad line objects), while red squares are LEG. When 
possible, we further mark LEG according to their FR type: 
crossed squares are FR II/LEG and dotted squares are FR I/LEG. 
The two solid lines represent the best linear fit obtained for the 
HEG and LEG sub-populations separately. The small black dots 
are the radio-sources from the SDSS/NVSS sample (see text for 
details) . The do t ted elh pse includes the Seyfert galaxies consid- 
ered bv lWhittlel(ll985h . 



the nearby source 3C 272. 1 (M 84) being fainter From the point 
of view of the narrow emission lines, the [O III] luminosity of 
BL Lacertae lies instead at the high end of 3CR/FR I. Thus this 
source is a strong outlier (by a factor of ~ 200-400 depending 
on the adopted value for its radio luminosity) from the relation- 
ship between line and radio luminosities followed by 3CR radio- 
galaxiefl This is true not only considering FR I, but also from a 
comparison of the LEG and HEG classes, and casts substantial 
doubts on the association of BL Lacertae with FR I, but more in 
general with the radio-galaxies in the 3CR sample. Note that also 
oth er classical samples of radio galaxies, like the B2 and the 2 
Jy ("Morganti et al.l ll997l:lTadhunter et"anil998h follow relations 
between radio and li ne emi ssion similar to the 3CR. 

iBaldi & Capeltil (|2009|) showed that a large ratio of line 
emission to the radio power with respect to 3CR sources is 
chara cteristic of the radio-loud AGN selected by iBest et aP 
( 12005 ). The latter authors cross-matched the ~ 212000 galax- 
ies drawn from the SDSS-DR2 with the NVSS and FIRSlQ ra- 
dio surveys, selecting a sample of 2215 radio-loud AGN (with 
a radio flux threshold of 5 mJy) to which we refer hereafter 
as the SDSS/NVSS sample. iBaldi & Capettil (1201 Ol) explored 
the spectro-photometric properties of the SDSS/NVSS objects 
showing that they are generally hosted by massive early-type 



* Other two strong outhers from the hne-radio correlation are present 
in the 3CR, namely 3C 371 and 3C 084; intriguingly, 3C 371 is a well- 
known BL Lac, while 3C 084 (although often classified di fferently) 
hosts a highly polarized and variable optical nucleus (M artin et alJ 
1197^ . typical of BL Lac objects. 

^ Sloan Digital Sky Survey, Data Release 2 dYork et alj I2OO0I) . 
National Radio Ast ronomy Observatory (NRAO) Very Large Array 
(VLA) Sky Survey jCondon et alJll99^ . and the Faint Imag es of the 
Radio Sky at Twenty centimeters survey dBecker et alj|199a) . respec- 
tively. 



galaxies with a low excitation emission line spectrum. From the 
point of view of its narrow line luminosity and extended radio 
power BL Lacertae falls well within the region covered by the 
SDSS/NVSS sample (see Fig. |5]l, although at a slightly higher 
L[0 III] than the region of higher galaxies density. 

However, these SDSS/NVSS sources do not show prominent 
broad lines like those observed in BL Lacertae. Because from the 
point of view of its radio flux "misaligned BL Lacertae" objects 
would be included in the SDSS/NVSS catalogue even if located 
up to a redshift of ~ 0.2, the only possibility before we con- 
clude that BL Lacertae is a unique object, at least in the nearby 
Universe, is that "misaligned BL Lacertae" objects were rejected 
on an optical basis. 

Indeed, iBest et afl (l2005h excluded from their sample 
sources recognised by the automated SDSS classification 
pipeline as QSO, because a bright nucleus prevents a detailed 
study of the host properties. These objects instead generally do 
show broad lines. Hence one may wonder whether it is possible 
that the parent population of BL Lacertae can be found among 
them. We have then selected from the DR7 of the SDS^I afl ob- 
jects classified as QSO with redshift z < 0.1, resulting in 1033 
sources, 90 of which have a radio source in the NVSS catalogue 
brighter than 5 mJy within 15". We then looked for objects with 
properties similar to BL Lacertae, filtering out 42 sources with r 
band and/or radio luminosities more than four times fainter than 
BL Lacertae, and 13 well-defined spiral galaxies, ending up with 
35 candidates. 

We downloaded their spectra from the SDSS public archive. 
Despite their automated classification as QSO, only 18 of them 
do show prominent broad emission linefl These galaxies, how- 
ever, have [O III] luminosity far higher than BL Lacertae, within 
the range L[0 III] ~ 2 x 10"*' - 4 x 10"*^ erg s"' (with a median 
of ~ 5 X 10"*' erg s"'), a factor of 5 to 100 higher than in BL 
Lacertae, and are also characterized by a HEG spectrum. These 
sources have line and radio lu minosities as well as line ratios 
typical of Seyfert galaxies (e.g.' Whittlel[T985rFI . 

Summarizing, we did not find any object in the SDSS/NVSS 
sample with isotropic properties similar to BL Lacertae. In par- 
ticular, only a few sources show broad emission lines and at the 
same time a luminosity of the host galaxy and a radio power at 
least a quarter of BL Lacertae; however, their narrow line lumi- 
nosities exceed that of BL Lacertae on average by a factor of 
12, and their HEG spectrum contrasts with the tentative identi- 
fication of BL Lacertae as LEG. Because, as mentioned above, 
"misaligned BL Lacertae" objects would be found in this ex- 
tensive catalogue, we must conclude that we failed to find the 
parent population of BL Lacertae, which apparently stands out 
as an "orphan" AGN. 



Note that the DR7 includes about four times more galaxies than the 
DR2 considered above. 

' For the remaining 17 galaxies, we set upper limits to their broad Ha 
luminosities fitting multiple emission line s to the spectra ob tained after 
subtraction of the stellar component, as in lButtiglione et aLl ( i2010h . All 
limits are between ~ lO**" and 10*' erg s"'. 

We remark one apparent exception, namely SDSS 
J122358.97+404409.3, with L[OIII] only 20% smaller than BL 
Lacertae, However, this object has a strong contamination by its active 
nucleus; the equivalent width (EW) of the NaD feature is EW(NaD) = 
1.04 A, strongly reduced by the dilution of its nuclear emission with 
respect to a typical value for early type galaxies of ~ 4 A, indicating a 
nuclear contribution of ~ 80 %. Its r band magnitude corrected for the 
nuclear contribution is Mf, = -20.7, nearly 2 mag below the host of BL 
Lacertae. 
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This result contrasts with the requirement that to each object 
observed with its jet (characterized by a buUc Lorentz factor F) 
forming a small angle with our line of sight a parent population 
of mis-aligned objects must correspond. Although in presence 
of nuclear obscuration not all of them would show broad lines, 
this population is expected to be formed by ~ p2 ^6 - 1000 
sources, having assumed a range of F ~ 4 - 30 from lJorstad et alJ 
(I2OOI . 

4. Summary and conclusions 

More t han a decade ago Verm eulen et alJ (fT99l and subse- 
quentlv lCorbett etal] ([1996. .2000 ) detected a broad Ha emission 
line in the spectra of BL Lacertae. This luminous line should 
have been detected in previous spectra, suggesting that its flux 
must have increased by at least a factor 5 since 1989. The lumi- 
nous Ha line suggested a Seyfert-like nucleus in BL Lacertae, 
complicating the already difficult task of understanding what 
AGN population this object (and BL Lac objects in general) be- 
longs to, considering that FR I generally do not show broad lines. 

To investigate this matter we acquired low- and high- 
resolution spectra of BL Lacertae with the TNG during four 
nights in 2007-2008, when the source optical brightness was 
R ~ 14-14.5. Our spectra confirm the presence of a luminous 
Ha broad fine of ~ 4 x lO"*' erg s"' and FWHM ~ 4600km s"', 
as well as several narrow emissio n lines. Through a cri tical com- 
parison of our data with those bv lCorbett et alj (l2000h . we con- 
cluded that the BLR luminosity has increased by ~ 50% in about 
ten years. This level of variability is not unusual for Broad Lined 
AGN and it does not necessarily implies an evolutionary trend. 

Then we examined the nuclear properties of BL Lacertae. 
The relationship between the SMBH mass and bulge luminosity 
in AGNs allowed us to derive a mass of 4-6 x 10^ Mq. Using 
the spectroscopic information to calculate the virial mass, we 
instead obtained a value about 20-30 times lower. To understand 
the reason of this discrepancy we analysed the disc and BLR 
properties of other AGNs, and concluded that the BLR of BL 
Lacertae is underluminous by a factor 70-300. 

Finally, we analysed the physical quantities that do not de- 
pend on orientation and beaming, and thus should also char- 
acterize the parent population of BL Lacertae. We defined di- 
agnostic indices with the most reliable narrow emission lines, 
and found that their values provide a tentative identification of 
BL Lacertae as a LEG. Broad lines are instead observed only 
in HEG, but the diffuse radio luminosity of BL Lacertae is at 
least 100 times lower than in these powerful radio sources. On 
the other hand miniature radio galaxies are LEG, share both 
the narrow line and radio power properties of BL Lacertae, but 
they do not show a BLR. Taking into account how the minia- 
ture radio galaxy sample was selected, we expect that it should 
include "misaligned BL Lacertae" objects, unless they were ex- 
cluded on the base of a QSO appearance. An analysis of the 
galaxy morphology, spectral features, and radio power of the 
QSO sources, initially discarded from the sample of miniature 
radio-galaxies, revealed that no object meets the requirements to 
represent the BL Lacertae parent population. Yet, for typical val- 
ues of the Lorentz factor, we would expect 10-10^ "misaligned 
BL Lacertae". 

This leaves us with the only possibility that the observed 
properties of BL Lacertae are the result of a transient short last- 
ing phase. We can envisage the following scenario, somewhat 
similar to that already suggested by Corbettetal. (1996). BL 
Lacertae in its initial state has properties similar to the sources 
of the SDSS/NVSS sample. Indeed, these are massive early-type 



galaxies and a large number of them have narrow lines and radio 
luminosities similar to that of BL Lacertae. From the point of 
view of their optical spectra they are LEG and lack broad lines. 
Subsequently (possibly ~ 20 years ago), its BLR underwent an 
increase of luminosity due to an increased amount of cold gas 
in the nuclear regions and/or to a higher level of ionizing con- 
tinuum. These two effects may even be related and caused by a 
fresh input of accreting gas. The BLR structure might not have 
yet reached a stable configuration, accounting for its different 
properties when compared to other AGN. Also the NLR lumi- 
nosity will grow with time and will also eventually change its 
state of ionization, but on a much larger timescale with respect 
to the BLR. 

Based on the analysis of a single object it is clearly impossi- 
ble to set a timescale for the duration of the putative bright phase. 
Furthermore, BL Lacertae was probably discovered since this 
object has been subject to repeated spectroscopic observations. 
However, our failure to find objects in the local Universe that 
might constitute its parent population suggests that the timescale 
associated with the period of high accretion must be orders of 
magnitude shorter than the lifetime of radio-loud AGN. 

An alternative possibility is that the birth of the BLR marks 
the transition from a low-power radio galaxy to a high-power 
source. This would require a rapid increase in the luminosity of 
the large-scale radio structures to reach the level observed in e.g. 
the HEG of the 3CR sample, within a sufficiently short time so 
as not to produce a substantial population of transient sources. 
Instead, the available data rule out that BL Lacertae became an 
AGN only very recently, i.e. that we are witnessing its birth, be- 
cause its radio emission extends ~ 10 kpc away from the core. 
This implies that this source is active since at least ~ 3 x 10^ 
years, assuming an expansion speed of 0.1 c. 

We conclude that the parent population of BL Lacertae can 
be found among the large population of miniature radio-loud 
AGN forming the SDSS/NVSS sample, but this also requires 
that this object is experiencing a short transient phase. A con- 
tinuation of the spectroscopic monitoring of this peculiar source 
caught in a crucial phase of its evolution can help us tremen- 
dously in our study of the physics and evolution of these sys- 
tems. 
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